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Abstract: 

The chemical compositions of stars from the Asymptotic Giant Branch are still poorly known due to 
the low temperatures of their atmospheres and therefore the presence of many molecular transitions 
hampering the analysis of atomic lines. One way to overcome this difficulty is by the study of lines in 
regions free from molecular contamination. We have chosen some of those regions to study the chemical 
abundance of the S-type star GZ Peg. Stellar parameters are derived from spectroscopic analysis and a 
metallicity of -0.77 dex is found. Chemical abundances of 8 elements are reported and an enhancement 
of s-process elements is inferred, typical to that of an S-type star. 
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1 Introduction 



The Asymptotic Giant Branch (AGB) is one of the 
most important phases of stellar evolution due to the 
rich nucleosynthesis that occurs during this stage, and 
therefore AGB stars are crucial objects for studies of 
Galactic chemical evolution. During this epoch stars 
produce important amounts of carbons and s-process 
elements which are brought to their atmospheres by 
deep dredged-up events and then ejected to the inter- 
stellar medium by stellar winds. Although fundamen- 
tal stones for the understanding of the final phases of 
low-mass stellar evolution and the stellar populations 
abundances, the chemical composition of such stars 
are still poorly known due to their low temperatures 
and to the blanketing of the atomic lines by molecu- 
lar transitions. Some authors have tried to overcome 
this problem by selecting spectral regions free from 
molecular contaminations or the so called molecular 
windows (e.g. Smith & Lambert 1985, 1990, Lambert 
et al. 1995, Vanture & Wallerstein 2002a,b, 2003), and 
the use of the near-infrared and infrared regions (e.g. 
Origlia & Rich 2005, Cunha & Smith 2006, Woolf & 
Wallerstein 2005a,b). Although working with a smaller 
number of lines, the chemical picture of these stars can 
still be delineated, provided that a careful inspection 
of the continuum is made and an appropriate model at- 
mosphere is adopted. The knowledge of the elemental 
abundances of such stars is an important tool for the 
understanding of the many steps followed by red giants 
along the M-MS-S-C spectral sequence, and their role 
in the chemical evolution of galaxies. 

In the present work we report the chemical abun- 
dances of GZ Peg, a S-type star. S stars show molec- 
ular transitions of ZrO and LaO, indicating an im- 
portant enrichement of s-process elements. Although 
rich in those elements, lines of the unstable s-element 
99 Tc (which has a half-life of ~ 10 6 years) are not ob- 
served in GZ Peg spectrum (Lebzelter & Hron 1999). 
GZ Peg is therefore classified as an extrinsic S-star, 
which is predicted to have inherited its s-process con- 



tent in a previous mass-tranfer event from an AGB 
compagnion. In intrinsic S giants, lines of Tc are ob- 
served and they are believed to have recently formed s- 
process elements and dredged-up to their atmospheres. 
The improvement of the quality of the stellar spectra 
and of the atomic constants and the availability of new 
model atmospheres allow now a robust analysis of the 
chemical patterns of such stars. This is the first of a 
series of papers about the detailed chemical abundance 
of s-process enriched stars. In this series the following 
questions are adressed: (i) What is the detailed chem- 
ical pattern of the sequence M-MS-S stars? (ii) Are 
there S stars with different s-process contamination 
as seen in Ba-stars? And if yes, what are the factors 
determining the intensity of s-process contamination?; 
(iii) Does the stellar population to which the AGB star 
belongs affects the s-process production? If yes, what 
are the effects of this relationship? 

2 Observations and Analysis 

Observations were carried out at the 1.52m telescope 
of ESO (European Southern Observatory), La Silla, 
in September 1999. The spectra were obtained using 
the FEROS spectrograph (Fiber-fed Extended Range 
Optical Spectrograph) with wavelength range 356 to 
920 nm and a resolution of R = 48,000. Reductions 
were performed using the DRS (online data reduc- 
tion system of FEROS). A subsequent reduction was 
performed using the CONTINUUM, RVIDLINES and 
DOPCOR tasks of the IRAF package. The mean sig- 
nal to noisy ratio of the spectrum is S/N = 61.97. 

Reddening has been estimated according to Chen 
et al. (1998) and bolometric correction has been in- 
terpolated from the grids of Lejeune et al. (1997). GZ 
Peg is a semirregular (SRA) variable according to the 
General Catalog of Variable stars (GCVS, Samus et al. 
2004). The light curve of the star at AAVSO (Ameri- 
can Association of Variable Star Observers) has been 
checked and a visual magnitude range of ± 0.4 is found 
for this star, therefore an uncertainty of 0.4 mag is ex- 
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pected. We have adopted the bolometric luminosity of 
M 6oi = -5.02 from Guandalini & Busso (2008). The 
temperature was inferred by the IRFM using the em- 
pirical calibrations of Alonso et al. (1999) and 2MASS 
colors. The resulting values are T(V-K) = 3502 K, 
T(J-H) = 3621 K and T(J-K) = 3454 K. These tem- 
peratures have been checked by the excitation equilib- 
rium of the Fe I lines (Fig. 1), and a best fit has been 
found for T e fi = 3600K. Although S stars are predicted 
to have masses in a range of 1.5 - 2 Mq (Jorissen et 
al. 1998), GZ Peg is an extrinsic S star, therefore a 
series of tests have been made to better constrain the 
mass of this star. Tests requiring simultaneously (i) 
the excitation equilibrium of the Fe I lines and (ii) the 
curve of growth fitting, were made with masses of 3 
Mq to 1.5 M© and the iteration procedure converged 
for a model with 2 M©. This value has been adopted 
to derive the trigonometric surface gravity also using 
the Hipparcos paralaxe. 

The stellar parameters have been inferred by re- 
quiring the ionization equilibrium of the Fe I and Fe II 
lines to derive the spectroscopic gravity and metallic- 
ity, and zero slope in the [Fe/H] vs. EW (equivalent 
width) plot to infer the microturbulence velocity. The 
resulting stellar parameters are the following: T e ff = 
3600 K, log g = 0.06, [Fe/H] = -0.77 dex and £ = 1.88 
kms -1 . The [Fe I/H] vs. EW relation is depicted in 
Fig. 2 and the curve of growth fitting is given in Fig. 3. 
The model atmospheres adopted here is a new version 
of the MARCS models atmospheres (Gustafsson et al. 
2008). We have selected 2 Mq spherical models and 
interpolated between grids to obtain the best result. 

A change of 50K in the temperature yields a change 
in gravity of 0.025 dex and a change of 0.4 mag in the 
bolometric magnitude yields a change in gravity of 0.16 
dex. We have determined the errors in stellar abun- 
dances due to stellar parameters uncertainties taking 
into account a change of A log g = +0.2 dex in gravity, 
A T cft = +50 K in temperature, A [Fe/H] = +0.15 dex 
in metallicity and A £ = 0.10 kms -1 in microturbu- 
lence velocity. The resulting errors are given in Table 
[T] The uncertainties have been added in quadrature 
and the total errors are given in the last column. 
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Figure 1: Excitation equilibrium of the Fe I lines 
for GZ Peg. 
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Figure 2: [Fe/H] vs. EW of the Fe I lines for GZ 
Peg. 



Table [5] we give the final abundances with the number 
of lines used in the analysis and the statistical errors 
(*/VN). 



Table 2: Chemical Abundances 



Element 


[X/H] 


S a 


Number of lines 


Fe 1 


-0.77 


0.04 


13 


Ti 1 


-0.42 


0.11 


9 


Ni 1 


-0.62 


0.13 


6 


Zr 1 


-0.19 


0.11 


11 


Cr 1 


-0.28 


0.32 


2 


Nd 2 


0.35 


0.13 


5 


Y 2 


-0.10 




1 


Ce 2 


0.56 


0.27 


2 


Ba 1 


-0.05 




1 



Statistical error a j\J~N 



Chemical abundances have been inferred by the 
analysis of the equivalent widths of the lines. The line 
list and the EW of the lines are given in Table [3] In 



3 Discussion 

A somewhat low metallicity of [Fe/H] = -0.77 has been 
found for GZ Peg. Abundances of [Ti/Fe] = +0.35 and 
[Ni/Fe] = +0.15 are derived, and a slightly high value 
for Cr, [Cr/Fe] = +0.48 is found (although this element 
is represented only by two lines with a high scater be- 
tween them). For the s-process elements, we have Zr 
I and Nd II abundances better defined (with a larger 
number of lines) than Y II, Ba I and Ce II. We have 
found s-process abundances typical of a S star for GZ 
Peg. In Figure [4] we plot our data with those of two 
S stars from Vanture & Wallerstein (2003), MS and 
S stars from Smith & Lambert (1986) and Ba-stars of 
Smiljanic et al. (2007) and Allen & Barbuy (2006). As 
can be seen from this figure, GZ Peg shows supersolar 
abundance ratios for [Ba/Fe], [Zr/Fe], [Y/Fe], [Nd/Fe] 
and [Ce/Fe], and an enhanced pattern for [Nd/Fe] and 
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Figure 3: Curve of growth fitting. 



[Ce/Fe] when compared to most of the s-enriched stars. 
The overabundance level is similar to that of the in- 
trinsic S star S Uma (Vanture & Wallerstein 2003), 
which has a metallicity of [Fe/H] = -0.3. 

In order to better understand the origin of the 
abundance pattern of GZ Peg, we have compared in 
Fig. [S]the resulting abundances with those predicted 
by Goriely & Mowlavi (2000) for an AGB star with 
metallicity [Fe/H] = -0.70 and mass of 2.5 Mq which 
has undergone 10, 30 and 50 dredge-ups. Surface abun- 
dances predicted for 10 dredge-ups are depicted by the 
solid line, for 30 dredge-ups by the dotted line, and for 
50 dredge-ups by the dashed line. Except for the Ba 
abundance, which has been derived from only one line, 
GZ Peg data are better described by the 30 dredge- 
up model. Such analysis may indicate that this star 
has probably accreted matter from an AGB star which 
has passed through an important number of dredge-up 
events. 



4 Conclusions 

A detailed analysis of the extrinsic S star GZ Peg has 
been performed. The spectroscopic stellar parameters 
derived for this star are: T e g = 3600 K, log g = 0.06, 
[Fe/H] = -0.77 dex and £ = 1.88 kms" 1 . The derived 
s-process abundances are in agreement with those of 
S-type stars or other stars with enhanced s-process 
abundances. Comparing the abundance level of GZ 
Peg with theoretical surface abundances of AGB stars 
from Goriely & Mowlavi (2000) we have found a good 
agreement with the 30 dredge-up model. Such anal- 
ysis indicates that the pristine AGB star which con- 
taminated GZ Peg has passed through a considerable 
number of dredge-up events. 
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Figure 5: [X/Fe] vs. Z for GZ Peg compared 
to theoretical surface abundances predicted for 
AGB stars of Goriely & Mowlavi 2000. Our data 
are compared to models with [Fe/H] = -0.70 and 
M=2.5M Q . The solid line represents surface abun- 
dances predicted for 10 dredge-ups, the dotted 
line for 30 dredge-ups, and the dashed line for 50 
dredge-ups. Error bars are the same as in Fig. U] 
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Table 1: Changes in chemical abundances due to stellar parameters uncertainties. In the 
last column we give the total error. 



Element 


A log g 


A T eff 


A [Fe/H] 


a <e 


Total error 




+0.2 dex 


+50 K 


+ 0.15 dex 


0.10 kins" 1 


S 


Fe 1 


0.09 


-0.03 


-0.01 


-0.05 


0.11 


Ti 1 


0.05 


0.04 


-0.03 


0.00 


0.10 


Ni 1 


0.07 


-0.05 


-0.03 


0.00 


0.09 


Cr 1 


0.04 


0.02 


-0.01 


0.00 


0.07 


Zr 1 


0.07 


0.07 


-0.02 


0.00 


0.11 


Nd 2 


0.08 


-0.02 


0.08 


0.00 


0.11 


Y 2 


0.09 


-0.02 


-0.05 


0.00 


0.10 


Cc 2 


0.09 


0.01 


0.05 


0.00 


0.10 


Ba 1 


0.02 


0.05 


0.01 


0.00 


0.05 
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Figure 4: Comparison of the chemical abundances of GZ Peg with MS, S and Ba stars. The reference of 
the symbols are the following: filled circle - GZ Peg, downward pointing triangles - S stars from Vanture 
& Wallerstein (2003), open stars - MS and S stars from Smith & Lambert (1986), squares - Ba-stars from 
Smiljanic et al. (20007). For the error bars we have adopted the highest error value between Tables [Hand 
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Table 3: Line List - Central wavelength of the lines, chemical element, excitation potential, oscilator 
strength, atomic data reference and equivalent width. 
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0.820 


-2.637 


01 


206.0 


8364.240 


Til 


0.840 


-1.756 


01 


319.2 


8377.860 


Til 


0.830 


-1.612 


01 


392.6 


8396.900 


Til 


0.810 


-1.779 


01 


380.0 


8412.360 


Til 


0.820 


-1.483 


01 


388.6 


8416.950 


Til 


2.240 


-1.034 


01 


133.7 


8417.470 


Til 


2.120 


-1.951 


01 


124.8 


8426.510 


Til 


0.830 


-1.253 


01 


468.0 


8438.920 


Til 


2.260 


-0.800 


01 


153.7 


8450.890 


Til 


2.250 


-0.903 


01 


106.8 


4385.240 


Nil 


2.740 


-2.037 


03 


97.3 


7414.500 


Nil 


1.986 


-2.549 


03 


182.7 


7422.300 


Nil 


3.635 


-0.251 


03 


133.7 


7525.140 


Nil 


3.635 


-0.546 


03 


91.6 


7555.600 


Nil 


3.847 


-0.046 


03 


107.4 



Atomic References: 
01 - Woolf & Wallerstein (2005b) 
02 - Va nture fc Walle rstein (2002b) 
03 - VALD (h ttp://ams.astro.univie.ac.at~vald71 ) 
04 - Vanture & Wallerstein (2003) 
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Table 3: Continuation 



Wavelength 


Element 




lo. 


cr gf 


Reference 


EW 


7574.043 


Nil 


3.833 


-0 


.533 


03 


104.3 


7437.138 


COl 


1.956 


-2 


.876 


03 


64.6 


7553.963 


COl 


3.952 


-0. 


.955 


03 


43.1 


7355.935 


CR1 


2.890 


-0. 


.290 


01 


208.0 


7400.226 


CR1 


2.900 


-0. 


.110 


01 


202.2 


7462.364 


CR1 


2.910 


-0. 


.040 


01 


213.6 


8348.283 


CR1 


2.710 


-1. 


,830 


01 


112.0 


8976.880 


CR1 


3.090 


-1. 


.110 


01 


156.2 


7439.890 


ZR1 


0.543 


-1. 


.810 


03 


92.9 


7479.532 


ZR1 


1.834 


-1. 


.380 


03 


24.8 


7551.493 


ZR1 


1.584 


-1. 


.360 


03 


36.9 


7554.780 


ZR1 


0.520 


-2. 


.280 


03 


111.6 


7562 129 

1 \j j-j . _L %j 


ZR1 


0.623 


_2. 


.710 


03 


46.9 


8070.115 


ZR1 


0.730 


-0 


.790 


01 


141.1 


8133.011 


ZR1 


0.690 


-1. 


.130 


01 


178.8 


8212.577 


ZR1 


0.650 


-1 


.320 


01 


167.2 


8305.987 


ZR1 


0.620 


-1 


.660 


01 


131.4 


8389.491 


ZR1 


0.600 


-1 


.760 


01 


148.7 


8584.281 


ZR1 


1.875 


-1 


.320 


03 


44.2 


7427.416 


ND2 


1.490 


-1 


.380 


03 


45.0 


7481.286 


ND2 


0.205 


-2. 


.900 


03 


30.8 


7577.496 


ND2 


0.205 


-2. 


.930 


03 


54.0 


8594.883 


ND2 


1.140 


-1. 


.860 


03 


33.9 


8643.478 


ND2 


1.200 


-1 


.750 


03 


42.3 


7947.597 


RBI 


0.000 


-0. 


.170 


01 


104.6 


7450.276 


Y2 


1.748 


-1 


.590 


03 


43.8 


8810.854 


CE1 


0.302 


-1. 


.665 


03 


18.4 


7486.557 


CE2 


2.548 


-0 


,394 


03 


23.5 


8772.135 


CE2 


0.357 


-2. 


,515 


03 


75.4 


7488.077 


BA1 


1.190 


-0 


.230 


04 


29.3 



